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ABSTRACT 

Context. If the gas in the evolving cosmic halos is dissipating energy (cooling) then due to the variation of the gravitational potential 

the dark matter halo also undergoes a compactification. This is well-known as Adiabatic contraction (AC). 

Aims. Complementary to the AC we investigate the resulting dynamical behaviour of the whole system if the backreaction of the AC 

of DM onto the gas is taken into account. 

Methods. In order to achieve sufficient high resolution also within the central halo region, we use a crude fluid approximation for the 

DM obeying the adiabatic contraction behaviour. Further, we restrict ourself to spherical symmetry and vanishing angular momentum 

of the studied matter configurations. The computations are done using a first-order Godunov type scheme. 

Results. Our results show that the dynamical interaction between gas and DM may lead to significant shorter collapse times. If the 

gas cools the dynamical behaviour of the whole system depends strongly on the shape of the initial density profile. 

Conclusions. Our findings indicate that for a certain mass range of halo configurations the dynamical interaction between gas and 

DM might be important for the halo evolution and must be taken into account. 
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1. Introduction 

During the past decade considerable progress has been achieved 
in understanding the structure formation process in the universe. 
In particular, the dynamics and evolution of the dissipationless 
component, i.e. the formation and evolution of the dark mat- 
ter (DM) halos was studied in great detail. The outcome of 
all the cosmological simulations is that the dark-matter compo- 
nent forms a nearly stable halo with an almost universal den- 
sity shape which is well approximated by the profile fit given by 
INavarro, Frenk & White (1 997)1 (NFW profile). 



However, the density profiles derived from the observed ro- 
tation curves of galaxies (in particular, from that of LSB galax- 
ies) indicate on much shallower profiles than the NFW profile 
or indicate even on a core structure (Marchesini et al. (2002), 
IGentile et al. (2004)[ |Kuzio et al. (2006)| |Salucci et al. (2007)| >. 
Therefore, the question arises whether the DM halo profiles have 
really an universal shape and in particular, which mechanisms 
or processes could then influence or change the DM density 
profile towards the observed ones. The structure of the inner 
DM density profile (cuspy or core like) is important with re- 
spect to various processes within the innermost parts of galax- 
ies and clusters of galaxies. The galaxy-galaxy lensing is a very 
promising observational method to map the dark matter distribu- 
tion in galaxy halos. For a successful reconstruction the knowl- 
edge of the density profile for a single DM halo is necessary, 
however ( Kleinhein rich et al. (2006)} . The innermost DM den- 
sity profile, i.e. how much concentrated the DM is within the 
central parts, is of particular interest if considering the conse- 
quences of possible dark matter annihilation and estimating the 
y-ray luminosity of Galactic halos (see, e.g., Stoehr et al. (2003)1 



Contrary to the DM, the baryonic fraction, (Mi, ~ 0.2 Mqm) 
is able to dissipate energy, may cool, and eventually forms a con- 
figuration much more concentrated than the surrounding DM. 
Due to the gravitational coupling of both components the DM 
responds dynamically on the compactification of the baryonic 
gas and the density shape of the DM is changing. That effect was 
first studied by Eggen, Lynden-Bell & Sandage (1962) using 
the model of adiabatic contraction (AC). Ze ldovich et al. (1 980) 
used the approach of adiabatic invariants to consider 
an analogous situation with respect to mass constraints 
on Leptons. A well elaborated form of the AC model 
was given by Blumenthal et al. (1986) 
later 



investigated and tested 
Cardone & Sereno (2005)] 



The model was 



by Ryden & Gunn (1987) 



Sellwood & McGaugh (2005) 



|Ascasibar et al. (2006)|[Diemand et al. (2007)) 



Vasiliev (2006) Dutton et al. (2007) and has been partly modi 
fied by Jesseit, Naab & Burkert (2002) In particular, it was con 
sidered within a cosmological context by[Gnedinet al. (2004) 
From the above studies, it is not clear however which processes 
are mainly influencing the final DM profile. |Gnedin et al. (2 004) 
have investigated how mergers and different cooling periods 
for the baryonic gas act on the DM dynamics. However, within 
mergers mainly the interaction of the DM particles is important. 
In the very inner region the steepest matter distribution defines 
the DM profile. It has been mentioned that the formation of 
the final DM profile might depend on the relation of cooling 
and dynamical times. The AC model assumes a relatively slow 
(adiabatic) change of the potential. Since the convergence of 
the simulations within r « rvmai is still n °t proven, it is 
difficult to make a final conclusion about the density evolution 
in the central region. However, this is the region which is most 
affected by any concentration of the gas in the result of cooling 
processes. 
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We want to consider the relatively artificial situation of a 
spherical symmetric matter distribution with fixed masses and 
limited spatial extension. The mass relation of gas to DM is fixed 
as 0.25. This ratio on average varies for the observed cosmic 
structures (galaxies, clusters, dwarf galaxies, etc.) very much 
since it depends on the particular conditions of the their forma- 
tion. It is also correlated with the total mass of objects. Thus, 
the assumed unique mass ratio of about 0.25 is taken in our con- 
siderations as an arbitrary value. In order to explore the dynam- 
ical interaction between the DM and gas we start from a sta- 
ble stationary configuration and then we let the gas cool. Thus, 
no effects of particular matter infall or merger-like mass growth 
is considered. To enhance the resolution considerably we use 
for the description of the DM a crude gas/fluid approximation 
assuming that the violent relaxation processes are fast enough 
to satisfy a nearly Maxwellian velocity distribution of the DM 
particles at each time and within each shell, provided a suffi- 
ciently coarse grained shell distribution is considered. A similar 
treatment of the DM was used before in |Chieze et al. (1997)| and 
Teyssier et al. (1997) Therefore, we describe at high accuracy 
the interaction of two gases/fluids which interact only gravita- 
tionally, and only the gas component is able to dissipate/cool. 
This allows to consider not only the dynamical response of the 
DM on the baryonic cooling as described by the AC models but 
vice versa the back-reaction on the gas dynamics, too. 

Recently, Conroy & Ostriker (2007)| have investigated which 
kind of energy input processes in a system comparable to clus- 
ters of galaxies may prevent the cooling catastrophe. Their out- 
come is that the common action of various heating processes 
may be able to stabilize the system over cosmic time. However, 
a fine tuning is necessary and the whole system is very sensitive 
to the parameters. They have considered the dynamics of the gas 
assuming the DM potential to be static. Our studies demonstrate 
that the regard of the gravitational interaction between gas and 
DM increases the tendency to an onset of the cooling catastro- 
phe. The paper is organized as follows: In the next section we 
introduce our basic assumptions and give the main equations for 
the considered system. The numerical treatment is given in sec- 
tion 3. Section 4 contains the results for our various considered 
cases where the time evolution of the DM halo+gas system is 
followed if the gas cools. A discussion and summary of these 
results are given in section 5. For better systematics, the details 
of the used numerical algorithm are described in the appendix. 

2. Basic equations 

This section reviews our main assumptions and the basic equa- 
tions describing the system of a dark matter halo with embedded 

gas. 



2.1. Dynamics of the baryonic gas 

The baryonic fraction of the halo is described as an ideal gas. The 
hydrodynamical equations for the density pc, the radial momen- 
tum density pc uc and the energy density pc £c = pc Uq/2+sg,t 
of a gravitationally bound gas sphere (we assume spherical sym- 
metry) are: 



d_ 
dt 

-K. (PCU G )+— — 
dt r l or 



o i a r 2 



1 d_ 
r* dr 

[r 2 p C u 2 c ] 



]-» 



dp G _ 50) 

-5— - ~Pc -5- 
or or 



j:(pcSc) + —T-[r 2 u c (p c s c +p c )\ = -p c u c - A (1) 



with the gravitational potential <t> given by the Poisson equation 
and the cooling function A. Pressure p, internal energy eq and 
density pc are related to each other by the polytropic equation 
of state: 



eg,t - 



Pc 



(y- 1)pg 



(2) 



The gravitational potential is determined by the density dis- 
tributions of the gas pc and of the DM pdm via the Poisson equa- 
tion: 



1 5 r 2 30 1 A ^\ 1 



(3) 



To consider typical energy dissipation in the gas we assume 
cooling by recombination as given in Blac k (1981)| We further 
assume that the gas consists of hydrogen, only. In that case we 
have 



A = 2,85- 10~ 27 7/ 1/2 

■ (5,914 - 0.5 In T + 0, 01 184 T 1/3 ) n e n H+ 



(4) 



where n e , n#+ denote the number densities of electrons and pro- 
tons and T being the temperature. We consider only collisional 
ionization and suppose ionization equilibrium 



^eh(T) n H n e = a H (T) n H + n e 



(5) 



where Teh(T) and ag(T) is the corresponding ionization and 
recombination rate, respectively (Black (1981)). 

Since the gas is assumed to consist of hydrogen only and 
and due to electric charge conservation the relation between the 
number densities n e , «#+ and n - pc/niH can be determined 
straightforward as: 



n e = hh+ = n 



^eh + a H 



(6) 



and (HJi can be written in terms of n and T only. 

In order to simplify the calculations, we restrict ourselves to 
the above cooling mechanism troughout. One should be aware 
however that within large mass systems like galaxies and clus- 
ters of galaxies the dominant cooling processes are helium and 
heavy element line cooling and Bremsstrahlung. In a realistic 
treatment, the cooling function is a composite of the various 
contributions the importance of which depends on the precise 
temperature and density conditions at each moment for the con- 
sidered volume. 



2.2. The dark matter dynamics - a fluid approach 

The dynamics of the dark matter is given by the collission- 
less Boltzmann-equation for the one-particle phase-space den- 
sity /(x, v, t) 



df 



df 



Z v '7j£ + Z 



d<& df_ 

dxj dvj 







(7) 



dt 



r 2 dr 



dr 



along with the Poisson equation for the mean gravitational po- 
tential O (see Eq. |3]l 

To obtain fluid-like equations of motion for the dark mat- 
ter we compute the moments with respect to the velocities 
vi up to the second order using the CBE as described in 
Binney and Tremaine (1988) and take into account that the ve- 
locities vanish at r — > 00. Our main approximation is to assume 
that the 3rd order moments vanish. In particular, this is true if the 
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velocity distribution is locally Maxwellian. Furthermore, we as- 
sume the velocity dispersion to be isotropic. Under these crude 
assumptions one gets the system of equations to be closed. In 
addition to the Jeans equation we get then an energy equation. 
If we assign formally the quantity 



_ l 2 

Pdm - -z Pdm °~dm 



PDM 0~DM,r 



(8) 



to a pressure then we have a full analogy to the gas 
equations. Making such an assumption we seem to sup- 
pose a kind of collision between the DM particles. However, 
if considering sufficiently coarse grained shell distributions 
for halos in n-body simulations then one observes nearly 
Maxwellian velocity distributions within each mass shell (see, 
e.g. |Hoeft, Mucket&Gottlober (2004)) . This is most proba- 
bly the result of sufficiently fast violent relaxation during 
the collapse of the considered systems. With other words, 
the mechanism of violent relaxation (Lyndon-Bell (1967), 
|Kull et al. (1997)} has to redistribute the velocities of the dark 
matter faster then the dynamical timescale. For the initial con- 
figuration, which we will start from, this process of relaxation 
leading to some equilibrium and to corresponding velocity distri- 
bution is assumed to be completed at the moment t — 0. Further 
on, the time-dependent evolution happens slowly enough that 
shell crossing is negligible. 

We want to use that approximation to study the interplay be- 
tween the gas and DM dynamics if energy dissipation in the gas 
happens. After transformation to spherical coordinates and re- 
stricting ourself to isotropic velocity dispersion (i.e., /3 - 0) we 
obtain the following system of equations: 



d_ 

at 



d~t PDM 



(Pdm udm) 
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(9) 



with the total energy density p D M sdm - Pdm vj:/2 + 3 pdm o~ 2 /2. 
Equations (0 together with the Poisson equation provide a 
closed system for the dynamics of the dark matter. As mentioned 
above these equations exhibit a strong analogy to the hydrody- 
namic equations as described before. Using the continuity equa- 
tion we can perform the energy equation in Eq.© to 



D . 

— log 
Dt & 



Pdm 



dm' 



d d 

— + udm-t- 
dt Or 



log 



PDM 



DM' 



= 



(10) 



The quantity t = PdmIo'tim represents the phase space density 



(see Taylor & Navarro (2001)) and the logarithm of r is a mea- 
sure for the entropy. Eq. (II Oil states that the entropy is constant 
along the flux lines, and therefore also r = const.. Thus, the 
assumptions on the velocity distribution of the DM yields an ef- 
fective equation of state (e.o.s.) which describes the process of 
adiabatic contraction of the DM sufficiently well. For the pres- 
sure analog p DM = Pdm<t dm we § et then the relation p D M <* Pdm 
formally identically to the adiabatic e.o.s. for an ideal gas. 

Note, for the static case the Eq. ( fTOb is identically fulfilled, 
and neither any restrictions on the e.o.s. nor any explicite expres- 
sion for the e.o.s. can be obtained. 



3. Numerical calculations 

3.1. The initial configuration of the DM-gas system 

For the initial distribution we chose a polytropic gas sphere (see 
Binney and Tremaine (1988)) of finite extension. We combine 
this with 

1. a spherical DM distribution obeying an effective e.o.s. with 
J dm = 5/3 

2. a dark-matter distribution following a NFW-profile 

According to the previous findings, for the dark matter and 
for the gas may be assumed 



Pdm - k-dm • P DM 



and 



Pg = k c - p' c 



(11) 



with the constants k D m an d k g . The initial profiles itself are com- 
puted for hydrostatic equilibrium and it's analog for the DM 

dpoM d<& dp c d$> 

— — = ~Pdm — respectively — — = -p G — - , (12) 
dr dr dr dr 

Together with the gravitational acceleration g - -d<S>/dr com- 
puted by Equation (O, one gets the following system of equa- 
tions 



T" Pdm 
dr 



3 Kdm 



2-f d 

PDM> J r PG 



7 k c 



2-y 

Pc 



dr 



(r's) 



■AnG (p D M +Pg) r 



(13) 



which can be easily solved by a Runge-Kutta algorithm. The 
central density of the gas and of the dark matter, the constants 
kdm and kq introduced above, and the ratio of specific heats for 
the gas y are free parameters. The central density gradients are 
set equal to zero. If we want the configuration of both, the dark 
matter and the gas, to be bound at the same finite radius, then we 
have to satisfy the condition 



k dm - Kg 



1G- 

Pg 



(»•=()) 



(14) 



This eliminates one free parameter. For a better interpretation we 
chose the halo mass M#, the cut-off radius of the profile Rh and 
the ratio of specific heats y instead of the above introduced free 
parameters . As an additional constrain we fix the ratio of total 
dark matter mass and total gas mass to be 0.25. 

As an alternative we consider the initial DM density distri- 
bution realizing a NFW density profile. In that case, the den- 
sity profile for the baryons is computed as above, but taking into 
account the given NFW-profile for the DM within the Poisson 
equation: 



Pnfw 



r(r + r s ) z 



(15) 



Here, the normalization for 6 is chosen in a way that the ratio of 
the total dark matter mass and the total gas mass is fixed as 0.25. 
The characteristic radius is set as an additional parameter. Thus, 
the equations for the the initial gas profile are: 



« g 2-7 

TPc = Pc 

dr y Kg 



dr 



{r 2 g) = -4nG (p N FW + Pc) r 2 



(16) 



In a second step pdm is determined to keep the system initially 
in hydrostatic equilibrium (Eq. (fT2l) ). 
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3.2. The dynamical evolution 

Once the system is forced to leave the hydrostatic equilibrium 
it's further evolution is determined by the full timedependent 
hydrodynamic and Poisson equations. In our case, the subse- 
quent cooling is responsible for this process. In order to compute 
the evolution of the system, we choose a Lagrangian scheme 
which is equivalent to the above Eulerian scheme. This has the 
benefit of conserving the mass automatically, and the treatment 
of the outer halo boundary is much simpler. We use the above 
discussed formal analogy to define a pressure for the DM as 

_ 2 

Pdm - Pdm o~ dm 

equations for the dark matter as for the gas. 



/3, This allows us for using formally the same 
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Fig. 1. Time fio when the central gas density reaches 10 times 
the initial central density as function of the initial boundary radii 
T -1 dr'. (17) at mass Mh - 1O U M . Note, for fixed mass the Rh is a measure 
for the mean density. 



The gravitation is the only coupling between the two kinds of 
halo matter, all other calculations may be performed indepen- 
dently for each component using identical routines. The intercell 
fluxes are computed using a first-order Godunov type scheme. 
We modified the original (second-order) MUSCL scheme pre- 
sented in | van Leer (1979)1 I n our case, due to the achieved high 
resolution the first-order scheme is sufficient, and encountering 
strong shocks is not expected. The Algorithm is in detail de- 
scribed in AppendixlAl 

The sources consisting of gravitational acceleration and 
cooling are calculated in each cell center. The gravitational ac- 
celeration at a given radius, is computed using the sum of the 
enclosed dark matter and gas mass: 



4jtG (m DM + m c ) 



The cooling function A t is calculated according to Equations 
© and ©. 

Instead of the physical quantities we use a system of dimen- 
sionless quantities based on the initial central gas density po and 
the velocity uq = Co/ yjyc , with co being the initial central speed 
of sound of the gas. The dimensionless quantities are computed 
as: p - plpo, u - u/uq, ... Radius and time are tranformed ac- 
cording to 



'0 



> /«g/(47tGp ) , to = r /u = Vl/(47rGpo) (19) 



The dynamical equations remain unchanged under this transfor- 
mation, and Eq.([T8l simplifies to: 

(m DU + m c ) 



Note, without cooling the system is free scalable and the actual 
dimensions of the halo have no impact on the dynamics. The 
cooling function has to be computed using the physical density 
and temperature and must then be transformed according to 



Ar 



roro 



At. 



4. Results 

We performed a series of computations starting with initial sta- 
tionary configurations as described above. First, it was proven 
that these configurations (if allowing for their time evolution) 
are stable against small perturbations as long no cooling takes 
place. Then we allowed the gas of the system to cool. This drives 
the system out of the initial hydrostatic equilibrium. The fur- 
ther thermal and dynamical evolution is governed by the system 
of Equations (Q]i, © and (|9). For large halo masses compara- 
ble with those of galaxies and more massive objects, the system 
always collapses within a period comparable with the free fall 
time of a sphere with mean density p mem ge ~ Mh/R 3 h - The Fig. 
(18) Elshows a measure of the collapse time in dependence on the ini- 
tial boundary radii at fixed gas + DM halo mass, i.e. as function 
of the mean density. 

If the contraction of the DM fraction is taken into account, 
the collapse time is marginally smaller. We want to get infor- 
mation about the action of the contracting DM onto the stabil- 
ity behaviour of the whole system. Therefore, we consider the 
behaviour of halos with masses of about M(, = 10 6 M o . Those 
masses are stable against a short cooling period leading to a fi- 
nite contraction whereas subsequent cooling leads also to an un- 
limited collapse. However, as can be seen below, in that case a 
period of quasi-stability or delayed contraction may occur. This 
behaviour is more sensible with respect to the dynamics of the 
DM. In that sense halos with w M(, represent a limiting case 
and the action of the co-collapsing DM can be studied in more 
distinctness. 

In order to study the influence of the initial density profile 
we have run the time-dependent computations for a set of dif- 

(20) ferent radii of the initial gas configuration. The NFW profile is 
formally extending till infinity but due to the r 3 dependence 
at the outer region it decreases very fast. We have located the 
border of the initial gas sphere well within the region where 
Pdm — * but varied it extension within a relatively narrow range 
(Rh.G - [20, 28, 34]pc). For the case of an initially "polytropic" 
DM distribution we can achieve identical outer boundaries for 
the DM and gas. In this case we have varied the initial radii as 

(21) R HG - [16, 20, 24]pc. The particular values for R^c were cho- 
sen to demonstrate the possible qualitatively different results. 
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NFW profile. (R„ = 20pc) 



for an initial NFW profile of the DM. 
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Figures [2] and [3] show the typical evolution of the density 
profiles for gas and DM for an initial NFW profile and a "poly- 
tropic" profile (with formally jom - 5/3) for the DM distri- 
bution. When the central gas density increases by 5 orders of 
magnitude the central DM density is increasing by one order of 
magnitude in both cases. The density profiles are steepening with 
time. 

It is obvious that the larger the initial radii are the smaller the 
initial central density is. The time evolution of radii containing 
a fixed amount of mass (r(0.5M),r(0.1M),r(0.01M))is given in 
the Figs. [5] and [6] for NFW and "polytropic" initial DM distribu- 
tions respectively. 

The time evolution of the central densities shows the most 
pronounced differences with respect to the chosen initial param- 
eters. Fig. |4] shows the time evolution of the central density, the 
central temperature and the absolute value of the cooling func- 
tion for the initial "polytropic" (labeled by A,B,C) and NFW (la- 
beled by D,E,F) dark matter density profiles. In each case three 
different initial gas density profiles determined by the radii of the 
halo boundary are considered (see corresponding labels at plot 
AandD). 



The initial boundary radii chosen as small as (Rh - I6pc 
and R H = 2Qpc) correspond to steep initial gas profiles. Due 
to cooling the collapse happens within a relatively short time 
interval. After an initial period of strong cooling, the cooling 
function A (see plots C and F) decreases due to the decline of 
temperature but rises again catastrophically if the central density 
increases. 

For shallower initial density profiles corresponding to Rh - 
2Qpc and Rh = 2%pc after a first short contraction phase, an 
intermediate delayed contraction phase occurs. In this case, the 
central density is lower when the cooling function rises again. 
The further increase in density is not sufficient to compensate 
the decline of temperature in order to keep the cooling function 
further increasing. Thus, the cooling function is declining before 
it rises again like above and exhibits some hump. 

In the case of even shallower profiles (Rh - 24pc and 
Rh - 34pc) the period of delayed contraction is significantly 
larger, (note the logarithmic time scale). Due to the initially even 
lower central density the cooling function exhibits a series of 
oscillations until a final unlimited raise of the cooling function 
occurs. 
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Fig. 4. Time evolution of central gas density (A,D), central temperature (B,E) and of the absolute value of the cooling function (C,F). 
The left column shows the quantities for the case of an initial "polytropic" DM distribution, the right column shows the behaviour 
for an initial NFW profile of DM distribution. 



In all cases, after a very short initial cooling period the tem- 
perature is varying only very little during the contraction phase. 

The oscillations in the plots for the shallowest profile appear 
as nearly periodic. Fig. [8] shows this in a blow-up of the plot 
during the delayed collapse for two different e.o.s. for the gas 
(y = [1 .4, 1 .666]). The quasi-periodicity depends on the value of 
y for the gas. 

For an interpretation of the above described behaviour we 
refer to Fig [7] This figure shows the evolution of the character- 
istic dynamic time t p , pressure response time t SO und ff l/ yT and 
cooling time t CO oi f° r the intermediate configuration (Rh - 20pc, 
initial "polytropic" DM configuration) in the central region. We 
notice that after a short cooling period the pressure is not longer 
able to balance the gravitational force: the dynamical time t p be- 
comes smaller than all other characteristic times and the system 



undergoes contraction without any pressure resistance. The gas 
falls almost freely. If the initial density is not too high the col- 
lapse gets delayed due to again increasing pressure force (the 
net acceleration becomes less than for pure gravitational infall). 
However, at nearly constant temperature the density enhance- 
ment due to contraction enforces the cooling process again (the 
cooling time decreases till a minimum) in a way that pressure 
resistance breaks down and the gravitational collaps will be en- 
forced. In the discussed case, this delay happens only once. Then 
the configuration is driven into the cooling catastrophe, i.e. both 
the cooling and the dynamical times approach zero. For even 
shallower initial density profiles, these periods of collaps delay 
may happen several times. 

We wanted to study the influence of any DM dynamics on the 
collapse behaviour of the gas. To this end we repeated the above 
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Fig. 8. A blow-up of the tome evolution of the central density Fig. 10. The same as in Fig. |9]but for an initial "polytropic" DM 
during the epoch of delayed collapse for two different y distribution. 



described calculations under the assumption that the distribution 
of the DM remains static, i.e. the DM appears as a source of 
an additional gravitational potential, only. This approximation is 
often used for only studying the dynamical evolution of the gas 
embedded in the DM halo (see, e.g., |Conroy & Ostriker (20 07)). 
In Figs. [9] we show the time evolution of the central den- 
sity in true units for initial NFW DM configurations with 
Rh = [20, 28, 34]pc (higher initial central density corresponds 
to smaller initial Rh). The upper panel shows the evolution of 
the DM central densities and the lower panel shows the evolu- 
tion of central gas densities. The dashed lines in the lower panel 
show the gas evolution for the case of static DM configuration. 
It can be clearly noticed that the evolution goes much faster, 
in general, if the dynamical evolution of the DM is taken into 
account. The shallower the initial density profile is the larger 
the amplification effect appears and the more extended the pe- 
riod of delayed collapse is. In Fig. [10] the analogous results 
for the case of an initial "polytropic" DM configuration is pre- 
sented (Rh - [20, 28, 34]pc). Note the appearance of the quasi- 
oscillations also for the evolution of the central DM density. 



A measure for the collapse strength could be the period of 
time fio during which the central density was increasing by a fac- 
tor of 10 with respect to the initial value po. The Fig. [TTI shows 
for different polytropic indices the dependence of fio on the cho- 
sen initial radii in the case of an initial NFW DM density profile. 
For given mass, the initial radii Rh are a measure for the steep- 
ness of the initial density profiles. The line style indicates on the 
chosen value for y, being the ratio of specific heats for the gas. 
Note, the effective yoM for the DM is fixed to be y^M - 5/3, 
throughout. 

It can be noticed that increasing Rh leads to an increasing 
collapse time. A larger y for the gas, i.e. a stiffer e.o.s. of the 
gas, leads to a larger t\o for given Rh- The figure shows pair- 
wise plotted graphs. In each case, the lefthand-sided graph shows 
the dependence in case of a static DM distribution. This means 
that taking into account a dynamical DM evolution coupled to 
the dynamics of the gas enhances the contraction, i.e. leads to 
shorter t\Q. One can distinguish three different segments in each 
graph. Each segment corresponds to a different type of collapse 
behaviour described above. The transitions between these differ- 
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Fig. 11. The collapse time fio is defined as the moment when the 
central gas density p reaches « 10p,„,. fio is shown as function of 
the initial halo radii (for a fixed mass Mh = 1O 6 M ). The initial 
density for the DM follows the NFW-profile. For each y the left- 
hand graph corresponds to a static DM distribution, while the 
right-hand graph corresponds to the dynamic DM evolution. 
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ent collapse phases is characterized by tight intervals of values 
R H . 

Fig.[T2lshows the dependences in the case of an initial "poly- 
tropic" DM density profile. Here, the static DM distribution 
leads to much longer t\o for given initial radii Rh in compari- 
son with the case of dynamical interaction between gas and DM 
(labeled by "Dynamic DM" in the plot). Interestingly, for the 
dynamic case the dependence on the polytropic index is inverse 
with respect to the above case of an initial NFW profile. Here, 
the three collapse phases seen above degenerate into two, but the 
transition between them happens within a much narrower range 
of Rh- In the case of the dynamic DM a break point at radii 
Rh ~ 22 - 23pc separates the short-time collapse from the de- 
layed collapse. 



5. Conclusions 

Our goal was to follow the time-dependent behaviour of a spher- 
ical halo consisting of DM and gas with high accuracy in a ID 
Model. In particular, we wanted to study the consequences of the 
dynamic interaction between the two kinds of matter onto the in- 
ner matter distribution if the gas is allowed to dissipate energy. 

We made the assumption that the DM is dynamically react- 
ing fast enough on any change of the overall gravitational poten- 
tial via processes of violent relaxation. This roughly justifies our 
crude fluid approximation for the DM which actually describes 
the adiabatic behaviour of a collisionless particle configuration 
under selfgravitation and being coupled to an embedded dissi- 
pating gas cloud via gravitational interaction only. 

Our results show that the dynamical interaction of gas 
and DM leads to an enhancement of the collapse be- 
haviour of the system (in agreement with the results of 
Kazantzidis et al. (2006)). The characteristic times for the con- 
traction process due to the dissipation (cooling) of the gas com- 
ponent are very sensitive with respect to the initial density pro- 
files although the configurations are all in a stable hydrostatic 
equilibrium initially. The steeper the central density profiles are 
the shorter the characteristic collapse times are. For halo masses 
of x M(, = 1O 6 M examined by us, two collapse phases for the 
gas fraction can be distinguished: After an initial short-period 
contraction follows a quasi-stationary state the duration of which 
depends on the initial density profile and possibly on the actual 
cooling function. 

We have considered the processes of collisional ionisa- 
tion and hydrogen recombination cooling, only. Therefore, the 
cooling function A vanishes effectively below T ~ 10 4 K. 
However, for halo masses of about M(, - 10 6 M o molecular 
hydrogen cooling may be important. Detailed simulations by 
|Yoshida et al. (2003)| and |Yoshida et al. (2006)| show that this 
leads indeed to further gas collapse and eventually to the forma- 
tion of first stars in the central region. In this context, the authors 
stress the point that the processes as cooling and cloud formation 
are mainly determined by the particular dynamics of the gravi- 
tational collapse. Thus, including molecular hydrogen cooling 
will lead to an even steeper gas density profile and hence is ex- 
pected to strengthen the collapse behaviour of the whole consid- 
ered configuration. 

Our result may also enforce the problem described by 
Conroy & Ostriker (2007)| Their outcome is that only a sensi- 
tive fine tuning of different cooling and heating processes is able 
to stabilize configurations like clusters over a time comparable 
with the Hubble time, at least. Taking into account the dynamical 
interaction between thermo-dynamically evolving gas and the 
DM halo such a fine tuning may be difficult to adjust. Especially 
within the central region of the considered configurations the dy- 
namic interaction between gas and DM may alter the situation 
considerably. 

Our computations do not indicate on any situation that the in- 
nermost DM distribution might get shallower due to interaction 
with the gas. Thus, including self-consistently the gravitational 
interaction between DM and gas is not yet sufficient to solve the 
discrepancy between the NFW density profile for DM halos in 
simulations and the indications on a core distribution from ob- 
served rotational curves of galaxies. The tendency of even steep- 
ening the inner DM profile by gas cooling may have influence 
on the investigations of DM annihilation signals. 

Whenever the cooling time approaches the dynamical time 
and both are becoming fast decreasing the cooling catastrophe 
cannot be avoided. Although we have considered only the exam- 
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Fig. 13. An example for the modified collapse behaviour if the 
cooling time strongly decreases with increasing optical depth. 
For comparison, the curve on the left hand shows the central 
density evolution for a optically thin medium. 



pie of cooling by recombination, this seems the outcome for all 
cooling functions <x p 2 (s. also Conroy & Ostriker (20 07)| l, i.e. if 
the density is increasing the cooling time always decreases. 

These considerations assume that any dissipated energy can 
leave the system immediately. If the medium gets sufficiently 
opaque with respect to interaction of the radiation with matter 
(e.g. photon-electron scattering) the radiation may need consid- 
erable time for carrying away the dissipated energy which may 
increase the cooling time considerably. To illustrate that we have 
introduced a factor oc exp(-« cr r A) in front of the cooling func- 
tion. Here A denotes some characteristic length, cr the relevant 
cross section and n the number density, e.g., of electrons. At 
large enough densities n the optical depth r = -ncrjA may ap- 
proach unity or may become even larger. In that case the cooling 
time is not longer decreasing at increasing density and one could 
expect that the epoch of delayed collapse may end at a station- 
ary configuration for particular cases. An example is given in 
FigH3]in comparison with the above considered optically thin 
case. This may indicate on a possible self-regulation process for 
halos within a certain mass range. Above all, this concerns dense 
objects as stars but possibly also the behaviour of the very cen- 
tral regions of forming cosmic structures. 

Our results demonstrate that the gravitational interaction be- 
tween gas and DM may have important influence on the dy- 
namical evolution of the halo systems during cosmic evolution. 
We have not included the possibility of continuous matter in- 
fall/accretion onto the evolving gas-DM halos. Furtheron, we 
did not include non-isotropic velocity dispersions and a non- 
vanishing angular momentum. Considering this may alter the 
situation (see, e.g., recent work by |Tonini et al. (2006)| l and will 
be subject to forthcoming research. 



Appendix A: Numerical Algorithms 

In the first step, velocity u* and pressure p* at the cell interfaces 
are determined by an exact Riemann solver. With these values 
the position of the cell interfaces after half of the timestep are 
computed: 



T i^ 



(t) + - At u* , 

2 '+5 



Because the mass coordinates of DM and gas are independent, 
the gravitational acceleration for the DM is computed by inter- 
polating the enclosed mass of the gas using the trapezoid rule: 
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with the gas cell j, being the outermost gas cell enclosed 
(\Tdm\i > WqIuL > l/DAfL'-i)' The procedure to calculate the 
gravitational acceleration for the gas is analogous. For computa- 
tion of the cooling function the temperature is computed accord- 
ing to 



T f = — Pi T; 



(A.4) 



Then, the cooling function [At]; can be calculated using 
Equations @, (O and (fJJJ. To conclude the timestep, the quan- 
tities at the time t + At are now computed (in this order): 
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The boundary condition at r — is realized by a reflecting 
boundary. Central pressure and velocity are computed using the 
solution of the Riemann-problem between the innermost cell and 
a ghost cell with the same density and pressure, but with a neg- 
ative velocity. The boundary between gas and vacuum is treated 
by solving the vacuum Riemann problem (see |Toro (1999)[ p. 
138). The solution consists only of one rarefaction wave end- 
ing at the contact discontinuity. Thus, the pressure is p* — 0., 
there. Using the Lagrangian version of the Generalized Riemann 
Invariant (see | van Leer (1979")] p. 105): 
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(u* - u) + (r*a* - Tfl) = 

y — 1 \ / 
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one gets for the velocity: 
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(A.6) 
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(A. 1) with the Lagrangian speed of sound a - Jypjr. 
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Instead of a sharp cut-off, the NFW-profile reaches zero only 
for r — > oo. We use a reflecting boundary condition there to in- 
hibit any artificial inflow or outflow. Since the outer boundary of 
the dark-matter profile is much more extended than for the gas, 
this choice should not affect the dynamics of the central region. 
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